Simulating merging binary black holes with nearly extremal spins 
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Astrophysically realistic black holes may have spins that are nearly extremal (i.e., close to 1 in 
dimensionless units). Numerical simulations of binary black holes are important tools both for 
calibrating analytical templates for gravitational-wave detection and for exploring the nonlinear 
dynamics of curved spacetime. However, all previous simulations of binary-black-hole inspiral, 
merger, and ringdown have been limited by an apparently insurmountable barrier: the merging 
holes' spins could not exceed 0.93, which is still a long way from the maximum possible value in 
terms of the physical effects of the spin. In this paper, we surpass this limit for the first time, opening 
the way to explore numerically the behavior of merging, nearly extremal black holes. Specifically, 
using an improved initial-data method suitable for binary black holes with nearly extremal spins, 
we simulate the inspiral (through 12.5 orbits), merger and ringdown of two equal-mass black holes 
with equal spins of magnitude 0.95 antialigned with the orbital angular momentum. 

PACS numbers: 04.25.dg, 04.30.-w 
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I. INTRODUCTION 

Although there is considerable uncertainty, it is pos- 
sible that astrophysical black holes exist with nearly 
extremal spins (i.e., in dimensionless units spins close 
to 1, the theoretical upper limit for a stationary black 
hole). Binary black hole (BBH) mergers in vacuum typ- 
ically lead to remnant holes with dimensionless spins 
X ~ 0.7 — 0.8 dHl], although if the merging holes are 
surrounded by matter the remnant's spin typically could 
be higher than x ~ 0-9 [H El- Black holes can reach 
higher spins via prolonged accretion [!, [f|: thin accre- 
tion disks (with magnetohydrodynamic effects neglected) 
lead to spins as large as x ~ 0.998 @, while thick-disk ac- 
cretion with magnetohydrodynamic effects included can 
lead to spins as large as x ~ 0.95 Even without ac- 

cretion, at very high mass ratios with spins aligned with 
the orbital angular momentum, binary black hole mergers 
can also lead to holes with nearly extremal spins [9h11|. 
There is observational evidence suggesting the existence 
of black holes with nearly extremal spins in quasars [l2| , 
and some efforts to infer the spin of the black hole in mi- 
croquasar GRS 1915+105 from its x-ray spectra suggest 
a spin larger than 0.98, though other analyses suggest 
the spin may be much lower p^HlSj . 

Merging BBHs — possibly with nearly extremal spins — 
are among the most promising sources of gravitational 
waves for current and future detectors. Numerical simu- 
lations of BBHs are important tools both for predicting 
the gravitational waves that detectors will observe and 
for exploring the behavior of nonlinear, highly dynam- 
ical spacetimes. Following Pretorius' breakthrough in 
2005 [la], several groups have successfully simulated the 
inspiral, merger, and ringdown of two coalescing black 
holes in a variety of initial configurations; however, all 
prior BBH simulations have been limited to spins of 0.93 
or less, which is quite far from extremal. The parameter 
X is a poor measure of how close a black hole is to ex- 
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FIG. 1. The rotational energy of a Kerr black hole as a 
function of the hole's dimensionless spin parameter \ := 
Spin/(Mass) 2 . The thick red line indicates the Bowen-York 
limit: standard Bowen-York puncture initial data — used in 
almost all numerical binary-black-hole calculations to date — 
cannot yield rotational energies more than 60% of the way to 
extremality. By using instead initial data based on two su- 
perposed Kerr-Schild holes ("SKS initial data"), in this paper 
we surpass the Bowen-York limit (green circle), opening the 
way for numerical studies of merging, nearly extremal black 
holes. 



tremality in terms of physical effects: a black hole with 
spin 0.93 has less than 60% of the rotational energy of 
an extremal hole with the same mass (Fig. [lj. 

Previous simulations have been unable to reach higher 
spins because of the way they construct their initial data. 
Just as initial data for Maxwell's equations must satisfy 
constraints (the electric and magnetic fields must have 
vanishing divergence in vacuum) , initial data for the Ein- 
stein equations must satisfy constraint equations. Most 



2 




Maom/M 0.9933 
7adm/M 2 0.6845 



d /M 15.366 
n M 0.014508 
do -0.0007139 



TABLE I. Properties of initial data evolved in this paper. 
The quantity M denotes the sum of the holes' Christodoulou 
masses at t = 0. Hole i (where i — A or B) has Christodoulou 
mass Mi and dimensionless spin \\ along the z axis (i.e., in 
the direction of the orbital angular momentum) . Also listed is 
the Arnowitt-Deser-Misner (ADM) mass Madm and angular 
momentum JI DM (e.g., Eqs. (25)-(26) of Ref. [H). The 
initial angular velocity fio, radial velocity do, and coordinate 
separation do were tuned to reduce the orbital eccentricity. 



BBH simulations begin with puncture initial data p7j . 
which assumes that the initial spatial metric is confor- 
mally flat (i.e. proportional to the metric of flat space). 
With this assumption, 3 of the 4 constraint equations can 
be solved analytically using the solutions of Bowen and 
York [IH, however, conformally flat initial data can- 
not describe single or binary black holes that are both in 
equilibrium and possess linear [2(| or angular [U,^ mo- 
mentum. Bowen- York puncture data can yield solutions 
of binary black holes with spins as large as x = 0.984 ini- 
tially, but when such initial data are evolved, the holes 
quickly relax to spins of about \ — 0-93 or l ess [HHHI- 
Several groups have evolved BBH puncture data with 
spins near but below this Bowen- York limit [|| l26l . [27ft . 
with the simulation by Dain, Lousto, and Zlochower [28[ 
coming the closest with spins of 0.967 at time t = 
quickly falling to 0.924. 

To reach spins beyond the Bowen- York limit, one must 
begin with initial data that is conformally curved. Re- 
cently, Liu and collaborators [2£| have constructed and 
evolved conformally curved initial data based on that of 
Brandt and Seidel [3(J HH for a single black hole with 
spins as high as x = 0.99. Hannam and collaborators [32 1 
have constructed and evolved conformally curved BBH 
initial data [HI, [34[ for head-on mergers of black holes 
with spins as large as x — 0-9. In Ref. [3]|, conformally 
curved BBH data with spins of x = 0.93 were constructed 
and evolved through the first 1.9 orbits of an inspiral, but 
no attempt was made to simulate the complete inspiral, 
merger, and ringdown. 

In this paper, we demonstrate that conformally curved 
initial data is suitable for simulations with nearly ex- 
tremal spins by using it to compute the first inspiral, 
merger, and ringdown of two black holes with spins larger 
than the Bowen- York limit. By surpassing this limit, 
our results open the way for numerical exploration of the 
gravitational waveforms and nonlinear dynamics of black 
holes that are nearly extremal. 



II. INITIAL DATA 

We evolve a low-eccentricity initial data set: a BBH 
where the holes have equal masses and equal spins of 



magnitude 0.95 antialigncd with the orbital angular mo- 
mentum. Some properties of the initial data used in this 
paper are listed in Table Q] 

Following Ref. 35| and the references therein, we 
construct constraint-satisfying initial data by solving 
the extended conformal thin sandwich equations with 
quasiequilibrium boundary conditions [36144 1| using a 
spectral elliptic solver [421] . The initial spatial metric is 
proportional to a weighted superposition of the metrics 
of two boosted, spinning Kerr-Schild black holes. 

We measure the quasilocal spin Sakv of each hole 
in the initial data using the approximate-Killing-vector 
method summarized in Appendix A of Ref. [35j . which 
is very similar to the prescription previously published 
by Cook and Whiting [43j . The dimensionless spin of 
each hole x is then related to Sakv by the formula 
X ■= S AKV /A/ C 2 hr , where M chr := y/ M^ + S 2 /4MT is 
the Christodoulou mass, M m := y/A/16n is the irre- 
ducible mass, and A is the area of the horizon. (For a 
single Kerr black hole, M c h r reduces to the usual Kerr 
mass parameter.) 

To reduce eccentricity, wc follow the iterative method 
of Ref. [44|, which is an improvement of the earlier 
method of Ref. [45|. For each iteration, we construct 
an initial data set and evolve it for approximately 3 or- 
bits. Then, the initial angular and radial motion of the 
holes are adjusted to minimize oscillations in the orbital 
frequency. Using this method, we reduce the orbital ec- 
centricity to approximately 10 . 



III. EVOLUTION 

We evolve our initial data using the Spectral Einstein 
Code SpEC . Building on the methods of Ref. H3 and 
the references therein, we have made several technical 
improvements to our code which both enable us to evolve 
our x = 0.95 initial data through merger and make our 
code more robust in general. Here wc briefly summarize 
some of the most important improvements; full details of 
these techniques will be described in a future paper. 

We use a computational domain with the singularities 
inside the horizons excised, and we use a time-dependent 
coordinate mapping to keep the excision boundaries in- 
side the individual apparent horizons as the horizons or- 
bit and slowly approach each other [48|. Our coordinate 
mapping also ensures that the excision surfaces' shapes 
conform to those of the horizons which enclose them. 
One important ingredient of our improved binary-black- 
hole evolutions is that the coordinate mapping is adjusted 
adaptively throughout the evolution, which is helpful be- 
cause the horizons' dynamics change from slow to fast 
during the simulation. 

Because we apply no boundary condition on the exci- 
sion surfaces, these surfaces must be pure-outflow bound- 
aries (i.e., must have no incoming characteristic fields) in 
order for the evolution to be well posed. A second im- 
provement to our code is that we now can adjust the 
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FIG. 2. Color online. Left panel: The z component \ z of the 
dimensionless quasilocal spin of one individual horizon vs time 
t. (The individual holes' spins are equal within numerical er- 
ror.) Right panel: spin of the common horizon vs. time and 
the final spins predicted by the fitting formulae in Ref. [4!|] 
("A"), Ref. Qj] ("B"), Ref. [U ("C"), and Ref. ("D"), and 
(for "B","C", and "D") the error bars corresponding to the 
fitting formulae's listed uncertainties. (Note that the horizon- 
tal positions of points "A"-"D" on the figure are arbitrary.) 
Our results are shown for several resolutions (labeled Nx or 
Nx.y, where x £ 3,4,5 and y G 1,2,3 label the resolution 
used before and after merger, respectively). 



velocity of each excision surface to keep the character- 
istic fields outgoing there. For the x — 0.95 simulation 
considered here, this characteristic speed control is nec- 
essary only during the last orbit before merger; earlier, 
it is sufficient to control the size of the excision surface 
using the method of Ref. [47j . 

A third element which we have recently added to SpEC 
is spectral adaptive mesh refinement. During the evo- 
lution, we monitor the truncation error of each evolved 
field, the resolution requirements of the apparent hori- 
zons, and the local magnitude of constraint violation; 
to maintain a desired accuracy, we then add or remove 
spectral basis functions as needed. In the simulation 
presented in this paper, we use spectral adaptive mesh 
refinement only during the final quarter orbit before 
merger. Throughout the entire simulation, we also adap- 
tively adjust the resolution of the apparent horizon finder 
as the horizon becomes more distorted. 



IV. RESULTS 

In Fig. [2j we plot the dimensionless quasilocal spin 
X measured on one individual horizon and also on the 
common horizon. From t = to t = 50M, there is a 
sharp, numerically resolved drop in the magnitude of the 
dimensionless spin x from 0.9498 to 0.9492. During the 
remainder of the inspiral, the spin drifts, with the amount 
of drift decreasing as resolution increases; at the highest 
resolution (iV5), the spin remains x = 0.949 throughout 
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FIG. 3. Color online. The orbital trajectory of the centers 
of the individual horizons and the individual and common 
horizons at the end of the inspiral (top left) and the real 
part of the I = 2, m = 2 mode of the emitted gravitational 
waveform h extracted at radius r — 405A/ (bottom). The 
holes travel through about 12.5 orbits before merging. All 
data is from resolution 7V5.3. 



the next 11.8 orbits, until just before merger, when the 
magnitude of the spin of each hole drops sharply. This 
result demonstrates that it is now possible to simulate 
BBHs where the holes retain spins beyond the Bowcn- 
York limit throughout the inspiral; it also opens the way 
for future explorations of the strong-field dynamics of 
merging, nearly extremal holes — dynamics that can only 
be explored using numerical simulations. 

During the ringdown, the spin x of the common hori- 
zon quickly relaxes to its final value of x = 0.3757±0.0002 
(where the uncertainty is estimated as the difference be- 
tween the highest and second-highest resolutions). This 
is approximately consistent with but slightly larger than 
the predictions obtained by extrapolating fitting formu 
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lae from simulations with lower initial spins of Ref. 
(Xfit « 0.371), Ref. Q fait = 0.372 ± 0.057), Ref. 
( Xflt = 0.369 ±0.012), and Ref. @ (xet = 0.366 ± 0.002). 
This result is a first step toward a better understanding 
of the relation between the properties of the remnant hole 
and those of the merging holes when the latter are nearly 
extremal. Numerical simulations that directly measure 
this relation (instead of extrapolating from lower-spin re- 
sults) will yield greater understanding of the properties 
of black holes produced from merging extremal holes. 

Figure [3] shows the individual horizon trajectories and 
the real part of the (t,m) = (2,2) spherical harmonic 
mode of the emitted gravitational waveform. We extract 
waves on a series of concentric spherical shells; the wave- 
form shown was extracted on the outermost spherical 
shell (at radius r = 405M). Accurate gravitational wave- 
forms obtained from this and future simulations with 
spins beyond the Bowcn-York limit will be useful for cal- 
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ibrating analytic template banks for gravitational-wave 
searches. 

The Christodoulou mass of the final black hole is 
M fina i/M relax = 0.9683 ± 0.0001, where M re i ax = 
1.0003M is the sum of the masses of the individual holes 
after the initial relaxation and where the uncertainties 
of Mgnai and M re i ax are estimated as the difference be- 
tween the highest and second-highest resolution. Under 
the assumption that each hole has a constant mass M ro i ax 
throughout the inspiral (which holds within O (l0~ 5 ) in 
our simulation after the holes have relaxed), the quan- 
tity 1 — Mfi na i /Mreiax represents the fraction of the initial 
mass that would have been radiated from t = — oo to 
t = +oo, had our simulation contained the entire inspiral 
instead of just the final 12.5 orbits. 

Our results demonstrate for the first time that it is 
possible to simulate merging black holes with spins larger 
than the Bowen-York limit of x = 0.93, the highest spin 
previously obtainable. Because astrophysical black holes 
may be nearly extremal, these simulations have astro- 
physical as well as physical relevance. In particular, this 
work opens the way to use numerical simulations to ex- 
plore the strong-field dynamics of merging, nearly ex- 
tremal black holes, to gain a better understanding of the 
properties of the remnant hole formed by a nearly ex- 
tremal BBH merger, and to provide high-spin gravita- 
tional waveforms for data analysis. 
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